Experimental phase-space-based optical amplification of scar modes 
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Waves billiard which are chaotic in the geometrical limit are known to support non-generic spa- 
tially localized modes called scar modes. The interaction of the scar modes with gain has been 
recently investigated in optics in micro-cavity lasers and vertically- cavity surface-emitting lasers. 
Exploiting the localization properties of scar modes in their wave analogous phase space represen- 
tation, we report experimental results of scar modes selection by gain in a doped D-shaped optical 
fiber. 



Wave chaos is investigating the wave behavior of 
bounded systems with chaotic classical dynamics. The 
common feature of chaotic modes in a linear (passive) 
medium is ergodicity [1]: modes of a chaotic cavity ex- 
hibit a generic statistically uniform field distribution Q . 
Scar modes that result from constructive wave interfer- 
ences along the least unstable classical periodic orbits 
break the ergodicity property [3] but are rather scarce 
events. Even though scar phenomenon has been inten- 
sively investigated since their numerical discovery by Mc 
Donald and Kaufman Q, their localization feature still 
arouses interest as in the profilic context of graphene in- 
vestigations 

Interaction between spatially localized scar modes and 
nonlinear effects is of fundamental interest and con- 
stitutes a recent evolution of Wave Chaos. The sub- 
tle interplay of gain with wave chaos and specifically 
scar modes has been widely examined in optics in 
micro-cavity lasers ®^ and in vertically-cavity surface- 
emetting lasers |{|. The problem of the interaction be- 
tween the lasing medium and the light field in a two- 
dimensional (2D) laser cavity is still a vivid research 
area [10]. In the presence of a saturable nonlinear effect, 
which does not involve any phase change, scar modes 
take a crucial part on the improvement of the lasing ef- 
ficiency. Thanks to their low-loss properties compared 
with generic ergodic modes, scar modes are commonly 
selected during the lasing process leading to enhanced 
emission directivity and lowered lasing threshold [ll| . In 
recent studies, the use of gain has also proved to favor 
the selection of scar modes in guided optics: numerical 
investigations revealed that the interaction of the chaotic 
modes with a localized active medium in a D-shaped op- 
tical fiber amplifier led to a selective amplification of scar 
modes [l2l ]. 

In this letter, we report the first, to our knowledge, ex- 
perimental demonstration of a selective amplification of 
scar modes through gain in a highly multimode chaotic 
system. This amplification is based upon the realiza- 
tion of a D-shaped optical fiber doped with a spatially 
localized gain. The amplified scar modes are built on 
the 2-bounce periodic orbit (2-PO) (see inset of FigJT]) of 
the associated D-shaped billard. Using the specific fea- 
tures of the scar modes in their wave- analogous phase 



space representation and a spatially coherent amplifica- 
tion process, we show that the selection mechanism of 
scars is strong enough to enhance scar modes and make 
them take over ergodic modes. Indeed, we can establish 
that certain non-scarred modes, displaying a good spatial 
overlap with the gain area, remain unamplified while the 
neighboring scar modes are clearly amplified, essentially 
because those non-scarred modes are not phase-space lo- 
calized in the way the 2-PO scar modes are. 

Our experimental system is a non standard Ytterbium 
doped multimode optical fiber with a core diameter of 
121 /mi, truncated at its half radius (see inset in FigJT]). 
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FIG. 1. Experimental setup for the amplification of scar 
modes in the double clad Ytterbium doped fiber amplifier. 
Li, lens (15- mm working distance), L2, lens (10- mm work- 
ing distance), L3, lens (2.5-cm focal length f), DM, dichroic 
mirror (T=90%@1064 nm, R=90%@978 nm), Pi, P 2 , 3-D po- 
sitioning supports, Fi, bandpass filter (3 nm @1064 nm), BC, 
beam collimator, BE, beam expander, M, mirror. Right inset, 
transverse cross section of the D-shaped fiber amplifier with 
the off-centered doped area (circle) and the 2-PO along the 
x-axis. The polymer cladding is not represented. Color online 

We fabricate the fiber preform by a standard MCVD 
(modified chemical vapor deposition) method. The core 
is made of pure silica including a 1 5- \im- diameter doped 
area of Ytterbium, free from any other co-dopant. A 
low refractive index polymer surrounds the silica core 
and acts as an optical cladding. The doped area is 
off-centered and located in a specific position deduced 
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from numerical investigations [12| which corresponds to 
the maximum of intensity for most of the 2-PO scar 
modes. The preform has been cut and shaped before 
being stretched for this purpose. The optical index mis- 
match between the doped area and the D-shaped core 
is induced by the ions concentration and takes the value 
An = 5 x 10- 4 ([Yb 3+ ] = 900 ppm). This value of 
An, coupled with a triangular index profile, prohibits 
any light guidance in the doped area. Both sides of the 
fiber are optically polished to avoid either at the input 
or at the output any diffraction induced by defects. 
The signal is given by a YAG laser (P s =400 mW) at 
A s =1064 nm and expanded thanks to a couple of micro- 
scope objectives (x4 and xlO) that play the role of a 
beam expander (BE). The resulting beam, of 5 mm of 
diameter, is sent through a 15-mm working distance lens 
(Li) and is transmitted by a dichroic mirror (DM, with 90 
% of transmission at A s ). The pump (P p =10 W) is pro- 
vided by a cw laser diode (coupled in an angle-polished 
fiber of 105-/im diameter) at A p =974 nm. The pump is 
focused thanks to a 10-mm working distance lens (L2) 
on the input of the D-shaped optical fiber through the 
dichroic mirror (90 % of reflexion at X p ). For an optimiza- 
tion of the pump absorption (excitation of a maximized 
number of modes), the beam is strongly focused into the 
multimode core [l3|, [l4| . 

The pump absorption coefficient has been measured un- 
der low pump power to prevent any population inversion 
and takes the value a p = 6.5 dB/m. The absorption and 
emission cross sections of the pump and the signal are 
respectively a pa = 2.65 x 10 -24 m 2 , a pe = 2.65 x 10 -24 
m 2 , a sa = 5.56 x 10" 26 m 2 and a se = 6.00 x 10" 25 m 2 . 
The gain of the signal has been measured along a 20 m- 
length fiber, and the maximum value G = 4.5 dB (with 
an output signal to noise ratio of 27 dB) is obtained at 
the optimal length L opt = 15 m. This gain value, which 
is rather low as compared to standards of telecom optical 
amplifiers, serves our purpose in demonstrating that am- 
plification is really at work in our highly multimode fiber. 
At the output of the fiber, the residual pump is filtered 
by a bandpass filter F\ of 3 nm-width around A s and a 
CCD camera is placed in the focal plane of a 2.5 cm-focal 
length lens to image the far-field intensity pattern. As 
we want to observe strong signatures of waves localiza- 
tion associated with scar modes, we choose to vizualize 
the signal intensity in the far-field configuration. Indeed, 
the square modulus of the spatial Fourier transform of 
the field distribution gives some information about the 
directions and modulus of the outgoing wave vectors. 

At the fiber entrance, we control the direction of the 
transverse wave vectors propagating in the fiber. A se- 
lection of the injected transverse wave vectors is obtained 
by a (f , ^-translation of the polished input of the opti- 
cal fiber in the divergent part of the focused signal beam. 
The figure [2^a) presents the typical far field intensity of 
the signal obtained at the output of the multimode 15 m- 
length D-shaped fiber amplifier when pump is off. Even 
if the injection of the signal is optimized to excite a scar 



mode among few ergodic modes, the bending of the fiber 
tends to distribute the signal intensity among numerous 
other modes. The observed statistically uniform distri- 
bution of intensity reveals the superposition of a large 
number of modes in the fiber. Figure EJb) represents a 
cut along the tv x -axis for K y = 0. It confirms the broad 
distribution of transverse wave numbers. 
For the same signal configuration, but with pump turned 
on, the figure [2fc) presents the far field intensity pattern 
of the outgoing signal. A pair of high intensity spots 
symmetrically distributed on both sides along the Ac^-axis 
appears in the center of the far field. These two peaks 
demonstrate that the signal propagates along the fiber 
with two privileged transverse directions given by k and 
— k along the x-direction. Those directions are associated 
with those of the 2-PO. The observation of two symetri- 
cal peaks along the 2-PO direction is the fingerprint of a 
scar mode in the far-field intensity representation as can 
be seen in the ft^-cut of the far-field intensity pattern of 
a typical calculated scar mode (Fig. [2ff)). 
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FIG. 2. Far field intensity observations (a), (c) and their 
transverse horizontal cut at k v — (b), (d) for the unpumped 
and pumped configuration respectively for k v — 7. Far field 
intensity pattern and transverse horizontal cut superimposed 
for the pumped configuration with k p = 6 (e) . Horizontal cut 
of the far field intensity of a calculated scar mode (f). Color 
online 

The transverse wavenumber corresponding to the ob- 
served experimental peaks can be estimated from fig- 
ure [2fd). The calibration of the far field representa- 
tion gives one pixel = (2.405±0.005) x 10 -3 /im _1 and 
then, the transverse wavenumber associated to the peaks 
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is deduced. We get k = (0.279 ± 0.003) /im" 1 = 
(17.0 ± 0.2) /R. This value has to be compared with the 
theoretical value deduced from a condition of phase co- 
herence of waves that bounce back and forth along the 
2-PO in the D-shaped billiard: 

n p C + Acj) + | = 2irp (1) 

where C = 3R is the 2-PO-length, Acj) is the phaseshift 
induced by the reflexion at the core/cladding interface, 
the additional tt/2 phaseshift is a consequence of the sin- 
gle self- focal point of the 2-PO, and p is an integer that 
gives the order of the scar mode. For p = 7, we get 
a value k p =7 = 17.08/i? which is in perfect agreement 
with the measured value: the scar mode of order 7 is 
selected thanks to the localized gain. We change the sig- 
nal launching condition by translating the fiber input in 
the (x, y) plane (Fig. [I]) in order to seed the amplification 
process on another scar mode. Figure[2fe) reports the en- 
hancement of two symetrical peaks along the ft^-direction 
for this new input illumination. The measured value of 
the transverse wave number ft = (15.0 ± 0.2) /R has to be 
compared to the calculated value (Eq. [1]) k p =q = 15.0 /R. 
The agreement is again remarkable and one can deduce 
that the scar mode of order p = 6 is enhanced thanks 
to amplification. These two examples are the proof that, 
by controlling the initial illumination condition, one can 
select a given scar mode along the 2-PO thanks to the 
optical amplification process. 

These preliminary experimental results demonstrate 
our ability to selectively amplify scar modes in a highly 
multimode fiber thanks to a localized gain. In agreement 
with common fiber amplification processes, the spatial 
overlap of the modes with the active area is usually the 
main control parameter of the amplification efficiency. 
We propose a complementary interpretation of the scar 
modes enhancement by gain based on the physical na- 
ture of scars. According to E. J. Heller, "scars manifest 
themselves as enhanced probability in phase space, as 
measured by overlap with coherent states placed on the 
periodic orbits" Therefore, scar modes, viewed as 
fingerprints of classical orbits in the near-field intensity 
pattern, also show strongly privileged wave directions as 
highlighted through their far-field intensity. Both spa- 
tial and directional features of scar modes can be most 
efficiently combined in the so-called Husimi representa- 
tion [15] . The Husimi representation is commonly used 
as a wave equivalent to the classical Poincare surface of 
section. Thus, it establishes a strong correspondence be- 
tween some particular modes and their associated clas- 
sical trajectories in the semiclassical regime. In the case 
of 2D cavities, the Husimi function measures the nor- 
mal derivative of the eigenfunction on the boundaries 
of the cavity [l6l-[l9|. Fig. [3fb-d-f) present the Husimi 
representation of three different calculated modes of the 
chaotic optical fiber with metallic boundary in the 2- 
dimensional space (#, ft sin x) with the curved abscissa, 
ft the transverse wavenumber of the modes and ft sin x 



the projection of the transverse wavevector along the 
tangential direction to the boundary (Fig. [3(a)). The 
Husimi representation (Fig. [3fb)) of a generic ergodic 
mode (Fig. [3(a)) generally explores a large part of the 
available space: no specific transverse wavevector direc- 
tions appear to be prevailing. This observation is con- 
sistent with the hypothesis of generic modes built upon 
an ergodic semiclassical behavior |2|, |20| . Ergodic modes 
result from the superposition of a large number of plane 
waves with fixed transverse wavenumber ft but random 
directions x- On the contrary, the Husimi representation 
(Fig. [3f d)) of a scar mode associated to the 2-PO (Fig. 
13(c)) exhibits a strong localization of the field around the 
angular abscissa 0o = and the exclusive direction of the 
2-PO, that is, xo = 0. All the 2-PO scar modes present 
such a similar localization attribute in their Husimi rep- 
resentation. 
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FIG. 3. Field intensity of specific modes and their respective 
Husimi representations. For a generic speckle mode (a),(b); 
for a scar mode along the 2-PO (c), (d); and for the mode 
with ftNS (e),(f). Color online 

The spatial and directional extensions of the Husimi pat- 
tern define a range of angular abscissa AO around 6q 
and of transverse wavevector directions Ax derived from 
A(ftsinx) around xo- Thus, a scar mode is built upon 
few highly directional plane waves, ensuring a strong spa- 



tial coherence. In a geometrical approach in terms of 
rays in the D-shaped billiard, this span of initial condi- 
tions extracted from the Husimi distribution determines 
an area in the billiard plane where rays are converging. 
This area is located in the vicinity of the so-called self- 
focal point of the 2-PO. This location of the active ions 
guarantees an efficient round-trip of rays through the ac- 
tive medium, as well as a spatially coherent amplification 
and also a good overlap with the field amplitude of the 
scar modes. We stress that all these conditions should be 
fulfilled to assure a scar mode amplification. Indeed, we 
observe that non-scarred (NS) modes can be insensitive 
to amplification despite a good spatial overlap with the 
gain medium. For instance, the field amplitude of the 
NS mode of Figj3^e) has a good spatial overlap with the 
doped area but is not selected during the amplification 
proces s 12111 . Using a Beam Propagation Method algo- 
rithm [22j, we simulate the propagation of a gaussian- 
like initial signal injected into the fiber with a transverse 
wavenumber k^s = 17.07 /R that corresponds to the spe- 
cific NS mode. The modulus of the transverse wavevector 
is fixed, but its direction is chosen so as to avoid the 2- 
PO direction and thus a 2-PO scar mode excitation. 
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FIG. 4. Intensity spectrum at the input of the fiber amplifier 
(dashed) and after amplification at the output of the fiber 
(gray) associated to a gaussian signal launched into the fiber 
with ftNS- Color online 

FigJH presents the intensity spectrum at the output of 
the D-shaped fiber amplifier. The targeted NS mode is 
not amplified whereas the scar mode of order p = 7 is 
strongly amplified despite its reduced contribution in the 
initial spectrum: the spatial field overlap is not able to 
clarify this favored amplification of the scar mode. The 
explanation is supplied in the analysis of the Husimi rep- 
resentation (Fig. EJf)). Indeed, the NS mode presents 
a broad distribution of transverse wavevectors directions 
in the Husimi representation that prohibits a coherent 
amplification mechanism. As a consequence, a strong 
localization in the Husimi representation around the sig- 
nature of the 2-PO is essential to ensure the efficiency of 
modes amplification. 

In this paper, we have presented an experimental ob- 
servation of scar modes amplification in the D-shaped 
optical fiber with localized gain. An active medium with 
a spatial location controlled by the phase space signa- 
ture of scars is a propitious system to selectively amplify 
scar modes. Investigations about scar modes in nonlinear 
systems are opening new attractive issues. For instance, 
having a scar amplifier may offer new perspectives for 
optical communications. The intensity spectrum of a 
scar modes amplifier has a response similar to the spec- 
tral response of a Fabry-Perot cavity. Each scar mode 
appears well-separated with a spacing between adjacent 
transverse wavenumbers equal to 2tt/C with kC! = 2^p. 
These well-resolved modes injected in a multimode fiber 
would be less sensitive to crosstalk and could be em- 
ployed as independent transmission channels in a process 
of mode division multiplexing p3j . 

We express grateful thanks to S. Trzesien and M. Ude 
for their priceless involvement in the preform and fiber 
manufacture. 
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